A Sinister Universe: Chiral gravitons lurking, and Lyth un-bound 
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Models of inflation involving non-Abelian gauge field backgrounds can produce gravitational waves 
at an observable level with a preferred handedness, while satisfying current observational bounds. 
This asymmetry comes about because the non-Abelian background generates parity-violation in the 
action for perturbations. In the specific model we study, Chromo-Natural Inflation, these gravita- 
tional waves can be produced at observable levels even when no held makes a super-Planckian fleld 
excursion, thus evading a common formulation of the Lyth bound. 



Scalar fields generically drive the inflationary epoch, 
Irowever some of the qualitative problems of inflation- 
ary model building can be ameliorated by incorporating 
gauge fields into the dynamics [IHS!- In particular, non- 
Abelian gauge fields can appear in a homogeneous and 
isotropic configuration, leading to slow-roll inflation even 
on steep potentials *! via competition between the po- 
tential forces and "magnetic drift" forces induced by a 
Chern-Simons interaction 'S]. In this Letter, we study the 
perturbations in such models, and show that the gauge 
field background results in a parity-violating action for 
fluctuations. During inflation, this parity violation im- 
prints on the spectrum of gravitational perturbations; 
one helicity is enhanced relative to the other, leading to 
a chiral spectrum of tensor fluctuations (chiral gravita- 
tional wave production was shown in a related context 
in El). 

The amplitudes of the gravitational waves (in ei- 
ther polarization) is not directly correlated to the (sub- 
Planckian) excursion of the inflaton field during inflation, 
even though this is a "single-clock" inflationary setting, 
because the motion of the inflaton along its potential is 
mainly resisted by "magnetic drift," not Hubble friction. 
Thus, this mechanism evades the simplest version of the 
so-called "Lyth bound" f^. A more comprehensive pre- 
sentation of our analysis will appear in a forthcoming 
work [F. 

Non-Abelian Gauge background: The Einstein- 
Yang-Mills equations are solved on FRW for a non- 
Abelian gauge field with this configuration: 

Ao^Q,A,^<i>5'iJa = a^5'lJa (1) 
Fo, - dr<i>5\r, = gcl>^ft^r. (2) 

where Ja is a generator of SU(2) [S]. This gauge field 
background violates parity by relating space and gauge 
indices with a particular orientation. Perturbing around 
tlris background, one finds parity-violating interactions. 
We focus mainly on the particular case of Chromo- 
Natural Inflation (hereafter CNI), with X an axion, in- 



teracting with the gauge fields via 

C=lR-l{dXr-V{X)-lTr[F^]-^XTr[FAF]. (3) 

In )1 we showed that this model inflates, with slow-roll 
provided via a magnetic-drift type force mediated by the 
Chern-Simons interaction. In the large drift force limit 
(A ^ 1), the slow-roll equations for this model are 

where g is the gauge coupling and / is the axion decay 
constant. Below we will take V = /L<^(1 + cos{X/ f)). 

Perturbations: Before delving into the details, let 
us sketch the novel features of the perturbations in this 
scenario. Because we have a non-Abelian gauge field on a 
background, we will have additional propagating scalar, 
vector, and tensor modes. The perturbative dynamics 
exhibits standard behavior at early times — fcr — >■ oo, 
where one recovers free plane waves; also, for late times 
kr — 0, the fluctuations freeze out, leaving only a single 
growing scalar and a single growing tensor mode. There 
are, however, two novel behaviors at intermediate times: 

1. For Am^ > '^ikT) > m^, where A — Xip/f, and 

= gip/H \s s. characteristic mass scale of gauge 
fluctuations, we find a complex interacting system 
for the scalars. The mode that governs the ampli- 
tude of the axion at superhorizon scales decays by 
a factor of order A. 

2. For a period kdr « around —{kr) « 1, one 
polarization of the gauge field tensor modes be- 
comes unstable and grows exponentially. This 
mode sources the corresponding gravitational ten- 
sor polarization, enhancing it relative to the other 
polarization. 

For the present work, we will ignore the vector perturba- 
tions, which do not play an important role. 

Some Details: We will work with the ADM metric, 

ds^ = -N^dT^ + a%jidx' + N'dt){dx^ + N^dt) (5) 
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where N is the lapse, A''* is the shift vector and hij is the 
metric on spatial hypersurfaces. We choose spatially flat 
gauge, in which (with jij transverse-traceless) 
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The gauge field perturbations decompose as 



(6) 
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(here t is symmetric and traceless) and the axion per- 
turbs to X + SX . In what follows, we work with explicit 
components of the fields by choosing the wavenumber 
along the direction. The gauge modes ^ then have 
an SVT decomposition in which = -^{tu ~ 122) ± iti2 
form a transverse traceless tensor, i3,(i±i2) a-nd Xi±i2 are 
transverse vectors, and z = ^(2^33 — tu — ^22), X3 s^ie 
scalars along with Sip. 

Tensors: The quadratic action for tensor modes is: 
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T2k9^{qj^ + tti^) + 2g^^{t^^^ + t^j^) (8) 
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where again 7± are the transverse traceless gravitational 
modes. The canonically normalized fields are thus 7* — 
7*/-\/2 and = y/2t^ . The crucial helicity dependence 
comes in the first term on the second line, which gener- 
ates a tachyonic mass over a range of scales near horizon 
crossing for t+ but not for i^, and leads to a substan- 
tial growth of the former near horizon crossing. Follow- 
ing this growth, the coupling between the gravitational 
and gauge modes (which also has a mild chirality depen- 
dence) feeds the disparate amplitudes of the gauge tensor 
modes into the gravitational sector. It is worth empha- 
sizing that this chiral splitting of tensor mode amplitudes 
is not solely due to the Chern-Simons interaction. The 
helicities are still split even in the absence of the Chern- 
Simons interaction (A — > 0), because of the preferred 
orientation implicit in the gauge field background. Any 
non-Abelian model relying on a non-trivial gauge texture 
of the form ([T]) will produce a chiral spectrum of tensor 
modes pO]. 

Working in the quasi de Sitter limit where H « const, 
a ~ —1/Ht, and (j) ~ Hcj), and introducing the variable 
X — —kr, the equation of motion for the gravitational 
wave can be written [TT] 
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X " Hx 
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where here and throughout primes denote derivatives 
w.r.t. X. The gauge field satisfies the approximate equa- 
tion 
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The equation of motion for t+ has an exact solution 

ikix) =AfcA/„,^(2zx) + BkWo,,0{2ix) (11) 

where Ma^0{2ix) and Wa^i3{2ix) are the Whittaker M 
and W functions, and we have defined 
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+ 2m^ = 2(1 + 2m|)/m^ = -2ia . (13) 



We set Ak and Bk, by demanding that the solutions 
approach canonically normalized positive frequency free 
plane waves as a; — > 00: 7*,£* — >■ e^^/V2k. We find 



A, = 



{2i)"C 
r(2/3 + l) 



r(. + /3 + i) ^''^ 



where C = r(— a + 13 + ^)/V2k. We can now solve for 
the resulting positive helicity gravitational wave mode 
by inverting Eqn. The resulting integrals can be 

evaluated in closed form in terms of Meijer G-functions. 
The late time result of this approximation scheme is 

7+ (a;) =7o(a;) + 2V'Bfc(/i -I- {h - TO^/3)) (15) 
where 70 (a^) is the homogeneous solution and 

(m^ — 2imrat + 2m — 2rrif^ sec (tt/?) sin (ttq;) T [a) 
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TT^ (jn^ + 2immt + 2m — 2TO( ) sec {n/3) esc {na) 
''2m{m + 2)r (a -I- 1) r (i - a - /3) r (i - a + /?) 
TT sec (tt/?) r (— a) 7r(m — iwt) sec (tt/?) 
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TT sec (tt/?) r (1 — a) 
mr(i - a-/3)r(i -a-f /3)' 
7r2(m -f imt)sec (tt/?) csc (ttq;) 



m(m + 2)r (a) r(i- a 
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We note that in the region where the approximation at 
(10) is breaking down (see the upper panel of Fig. [I]), 
the metric fiuctuations have frozen out and as a result, 
this approximation is excellent [I2j . 

Scalar perturbations: We must also give a brief ac- 
count of how scalar perturbations behave in this theory. 
A suitable gauge choice is non-Abelian Coulomb gauge 



A*' =5,**-i5[^„r] = 



(17) 



30 
20 
10 


-10 
-20 
-30 



1000 



100 



\ 

I \ 

— \ 

I \ 

- 4r]/H ^ 




S[r]/H ~^ \J\j 




/ 

: / 




/ 




/ 

/ 

/ 

~- i 1 1 < < 


\ / 5e[?+]/2 : 

iiV 1 , ; 


0.001 0.01 0.1 
X — — kr 

-1 


1 10 




FIG. 1. The upper panel, dashed hnes show the numerical in- 
tegration of the equations resulting from the action, Eqn. (Isl. 
Solid lines show the result of the approximation of Eqns. 
and The values of the parameters are chosen so that 

~ 2.1 and ip « 0.048. The lower panel shows the ratio of 
the power in the two gravitational wave helicities. 



which results in a Gauss law constraint that is algebraic 
in fc-space. This gauge condition can be used to eliminate 
X3 and, combined with the solution to the Gauss law 
constraint, leads to an action of the form 



is then determined in terms of the gauge scalars as 



X 



(19) 

Substituting this solution back into the action and mak- 
ing the further substitutions 

x(2m,^ Lp - X z) Am^xcp + {61711 + x'^)z 



3(2to2 



6(2mj 



(20) 



we find a simple reduced action depending only on 
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Note that the <p dynamics is unstable when < 2. 
In contrast to the tensor gauge modes, this scalar mode 
instability is present for a long time ^ x ^ and 
fully destabilizes the dynamics. There is thus a narrow 
range of that could be consistent with data: small 
leads to unstable scalar modes, while large leads 
to unacceptably large tensor modes. 

The amplitude of the axion on super-horizon scales is 
determined (up to order one corrections) by the behavior 
of if near a; = 0: 
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Gij^i<^[, + Mij<S>i^j + y^/,7$/$j (18) 



in terms of fields $ = {X, &/), z), with X = a5X the canon- 
ically normalized axion fluctuation; here the coefficient 
functions G, -M and A are rational functions of x with 
coefficients that depend on and A. The action is 
that of a charged particle moving in three dimensions 
in a time-dependent harmonic well and magnetic field. 
The important terms in the action at large A and x are 
a potential term X^m'^X'^ /x'^ and "magnetic" coupling 
(Xm^/x)(z — &j))X' . In the asymptotic past x ^ oo, the 
mass matrix — >■ 1, and .4 — ?> 0, yielding free field 
equations for the modes. For x < Am^, the mode of in- 
terest is a 'magnetic drift' mode where (as in the solution 
for the background) the large potential term balances 
against the large "magnetic" coupling, and the X kinetic 
term is largely irrelevant. A simplification of the dynam- 
ics ensues if we drop the X kinetic term altogether, as 
well as the subleading 0(A°) terms in Mxx- The axion 



The limiting value y(0) is of order one for large m^, 
starts growing around ~ and diverges as — 
due to the softening potential for this mode. 
We can now estimate the amplitude of scalar and ten- 
sor fluctuations. In spatially flat gauge, the curvature 
perturbation is 7?, = HSu, where 6u is the perturbation 
to the velocity potential, determined from 

To. ~ iPx + PYM)a^N, + aXd,5X (22) 
- (i/M(2z + 45/.)-#3_a^5Y) -252^Ar, . 

The term oc gip^X/f is much larger than the rest. It 
conies from the gauge field momentum 



To, = Tr(£; X B), « g^TTi[Aj,'^o][Aj, A,]) = g^^s^* 
where ^'g is determined by the Gauss law constraint 

A 
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(23) 
(24) 
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At late times, the leading piece of the first term comes 
from g'^[Ai, [Ai, V^o]], and the leading piece of the second 
term is {gX/ f){diSX)[Aj , Ak]e^^'' ; thus one finds 
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leading to the relevant term in (22 1. The scalar curvature 
perturbation that results is 
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where en = —H/H^ and an overdot denotes a deriva- 
tive with respect to cosmic time. The background has 
~ Ve/i^gXH) and H ~ VgX/{6fH), where Ve is 
the derivative of V w.r.t. 6 = X/f, thus 



5X 
X/H 
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Thus the axion is acting as the clock, as anticipated in 
[T], and the perturbation analysis is simply finding the 
magnitude of this perturbation. Now, relative to usual 
inflationary result for scalars on de Sitter, the magnitude 
of the SX perturbation is reduced from by a factor of 
A; also, the axion velocity is reduced by a compensating 
factor due to "magnetic drift". 

Putting it all together, one finds the scalar curvature 
fluctuation is given by 



V2{l + ml)\ Ve ) 



(28) 



(note that in the reduced system, the scalar d.o.f.s are 
normalized to unit amplitude in the far past); expressed 
in terms of the slow roll parameters 
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Relative to the corresponding expressions for single-field 



inflation, the slow- roll parameters (29) directly exhibit 



the additional suppression of motion along the inflaton 
potential by the coupling A to the gauge sector; and there 
is an additional enhancement from cp(0) when is small 
enough. Differentiating (30 1 with respect to dN ~ —Hdt, 



one finds the spectral index 



-lc:i-2eH + VH + d\ogcp{0)/dN . (31) 



The tensor-to-scalar ratio for the unenhanced graviton is 



{l + ml)eH 



"4, 



(32) 



while the other tensor-to-scalar ratio r+ is enhanced by 
the effect discussed above, determined from eq. (15 1. 



These formulae demonstrate that it is possible to satisfy 
current observational bounds with this model, while pro- 
ducing observable gravitational waves even in the absence 
of the chiral enhancement of one helicity mode. The chi- 
ral enhancement will be observable but not catastrophic 
for < 2 (sec Fig. [I]). 
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Note Added: While we were writing up our results, 
parallel work appeared |13j independently deriving the 
stability condition > 2 using a WKB method. 
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